WHAT IS CLAIMED IS: 
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1 . A method of removing streaks from multi-band digital images in 
an original spectral space, said method comprising the steps of: 

a) selecting a plurality of bands of a multi-band image to perform the 
streak removal operation upon; 

b) transforming each of the bands of the multi-band image to a spectral 
space advantageous for streak removal for multi-band imagery; 

c) performing a streak removal operation on each band in the 
advantageous spectral space using information from the other spectral bands; and 

d) transforming the streak removed bands from the advantageous spectral 
space back to the original display space. 

2. The method as claimed in claim 1 wherein the spectral space 
advantageous for streak removal for multi-band imagery is dependent on at least 
one of the number of the bands of data, the specfral bandpass of each of the 
imaging bands, and imaging band dependent characteristics of the one or more 
sensors used to capture the bands. 

3. The method as claimed in claim 1 wherein the transformation 
performed in step b) is a linear combination of the original bands. 

4. In a method of removing columnar streaks from a multi-band 
digital image of the type in which it is assumed that pixels in a predetermined 
spatial and spectral region near a given pixel are sfrongly related to each other and 
employing gain and offset values to compute streak removal information, the 
improvement comprising: transforming the pixels of the multi-band image to a 
spectral space advantageous for streak removal, wherein the transformation is a 
Unear combination of at least two of the original bands, and testing for a strong 
relation between the pixels in a predetermined spatial and spectral region near a 
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given pixel and computing streak removal information only if such a strong 
relationship exists, whereby image content that does not extend the full length of 
the image in the column direction will not be interpreted as a streak. 

5. A method of removing streaking in a multi-band digital image in 
an original spectral space, said method comprising the steps of: 

a) transforming the multi-band image to a spectral space advantageous for 
streak removal for multi-band imagery, thereby forming a transformed image; 

b) detecting pixel locations in the transformed image where pixel-to-pixel 
differences caused by streaking can be distinguished from normal variations in the 
scene data; 

c) performing a linear regression to determine an initial estimate of the 
gain and offset values between each pair of adjacent pixels in a direction 
perpendicular to the streaking using the pixel values at the detected locations; 

d) performing a statistical outlier analysis to remove the pixel values that 
are not from streaking; 

e) performing a linear regression to determine the gain and offset values 
between each pair of adjacent pixels in a direction perpendicular to the streaking 
using the pixel values at the detected locations that are not statistical outhers; 

f) setting the slope value to unity if it is not statistically different from 

unity; 

g) setting the offset value to zero if it is not statistically different from 

zero; 

h) using the slope and offset values to remove streaking from the 
corresponding line of image data; and 

i) transforming the streak removed fransformed image from the 
advantageous spectral space back to the original display space. 

6. A method for removing columnar sfreaks in a multi-band digital 
image, comprising the steps of: 
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a) transforming the multi-band data to a spectrally advantageous space, 
thereby forming a transformed digital image; 

b) selecting first and second adjacent columns of pixels from the 
transformed digital image; 

c) forming a colunm of pixel value pairs, representing the pixel values of 
the adjacent pixels in the two columns; 

d) forming colunms of local mean values, representing the mean values of 
pixels in an N-pixel window for each column; 

e) forming columns of mean-reduced values, representing the pixel value 
minus the corresponding local mean values in each colunm; 

f) calculating the correlation between the bands in the local region; 

g) forming a column of difference metric values, representing the sum of 
the squares of the difference between corresponding mean reduced values in an N- 
pixel window; 

h) forming a first reduced column of pixel value pairs by removing from 
the column of pixel value pairs, those pixel values whose absolute difference 
between the pairs is greater than a predetermined difference threshold; 

i) forming a second reduced column of pixel value pairs by removing 
from the first reduced column of pixel value pairs, those pixel values whose 
corresponding difference metric values are greater than a predetermined difference 
metric threshold; 

j) forming first slope, offset, and standard error values by performing a 
linear regression between the pair of pixel values in the second reduced column of 
pixel value pairs; 

k) forming a column of hnear prediction values using the slope and offset 
values and the first pixel value of the pair of pixel values in the second reduced 
column of pixel value pairs; 

1) forming a column of regression error values, representing the difference 
between the second pixel value of the pair of pixel values in the second reduced 
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column of pixel value pairs and the corresponding linear regression prediction 
value; 

m) forming a third reduced column of pixel value pairs by removing from 
the first reduced column of pixel value pairs, those pixel values whose 
corresponding regression error values are greater than a predetermined regression 
error threshold related to the standard error value; 

n) forming second slope and offset values by performing a Hnear 
regression between the pair of pixel values in the third reduced column of pixel 
value pairs; 

o) setting the second slope value equal to unity if it is determined to not be 
statistically different from unity; 

p) setting the second offset value equal to zero if it is determined to not be 
statistically different from zero; 

q) adjusting the value of each pixel in the second column of pixels in the 
digital image by multiplying each value by the second slope value and then 
subtracting the second offset value; 

r) repeating steps a-o for all adjacent columns of pixel values in the image; 

and 

s) transforming the multi-band data back to the original specfral space. 



